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From photometric observations of elliptical galaxies, among which are both radio galaxies and 
radio-quiet objects, an investigation was carried out of the relationship 'redshift - age of the 
stellar system' (Az/At). By means of this relationship cosmological parameters H{z) and f^A 
are estimated. Ages of stellar systems are determined within the framework of evolution models 
of synthetic spectra PEGASE and GISSEL. This approach can be considered as time study of 
objects of the early Universe independent of other cosmological models. Construction of a pooled 
sample is described, containing 220 objects from different populations of elliptical galaxies, for 
which an analysis of the upper limit of the age of formation of a stellar system was performed. 
These data were used to estimate the boundaries of determination of the cosmological parameters 
H„ and A-term: Ho = 72± 10 and Qa = 0.8 ± 0.1 in the model GISSEL and i/o = 53 ± 10, and 
ilA = 0.8 ± 0.1 in the model PEGASE. 
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1. Introduction 

An essential point in determining cosmological pa- 
rameters is the independence of the procedure of the 
classical methods, which have come to be classical, 
such as deep three-dimensional surveys of galaxies 
and clusters, la-type supernovae (see, for instance, 
Leibundgut 2001) or relic radiation (for instance Ef- 
stathiou et al., 2002). 

One of the independent techniques is based on the 
datings associated with the age of galaxies (for in- 
stance, Saini et al., 2000). The first attempts to esti- 
mate 51a with the use of ages stellar systems were 
made a few years ago (see, for instance, Parijskij, 
2001) proposed by Jimenez and Loeb (2002). It is 
based on the datings connected with variations of 
ages of galaxies determined by the spectroscopic tech- 
nique. 

Such an approach makes it possible to construct 
an independent chronological scale applicable to the 
early stages of evolution of the Universe. This ap- 
proach is based on measurements of differences of 
the ages Deltat between two passively evolving galax- 
ies which form at the same time but separated by a 
small interval Deltaz. Then one can determine the 
finite difference {Az/ At) k, dz/dt. All the galaxies in 
the procedure proposed by Jimenez and Loeb must 



have alike metallicities and low rates of star forma- 
tion (i.e. red colors), while the mean age of the sys- 
tem must be considerably larger than the difference 
of ages of galaxies At. By applying this differential 
method, Jimenez and Loeb (2002) suggest H(z) and 
LOQ (z) to be measured directly from the first and sec- 
ond derivatives (Az/At) {A^z/At^): 
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By proposing this differential method, the authors 
show that it is necessary both to increase the sample 
of galaxies and to improve the signal/noise ratio. 

We have used a similar method but with another 
type of dating of the age, namely, from photometric 
data and choosing optimum consistency of the distri- 
bution of energy in the spectrum (SED), depending 
on age, by the observed fluxes. This procedure, which 
has already become standard, operates with sufficient 
stability for a pure sample of elliptical galaxies (see, 
for instance, Verkhodanov et al., 1999), though it may 
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give an error in the age up to 2 Gyr. 

The methods using color and spectral ages of 
galaxies arc based on chronometrying the rates of 
expansion of the Universe from physical processes 
not associated with cosmology: from the rates of nu- 
clear reactions in stars, the knowledge of which for 
standard stars of the solar type are accurate enough 
and in the last decade have obtained numerous di- 
rect and indirect confirmations, including the latest 
achievements of acoustic tomography of the entrails 
of the Sun. For this reason, the proposed method of 
chronometrying the evolution of the Universe resem- 
bles those of chronometrying on Earth from the data 
of radioactive decay and in any case is independent 
of others. 

Note that giant elliptical galaxies with high radio 
luminosity and with the old stellar population are 
the most suitable objects for estimating the age of 
stellar systems. The present-day models predict fast 
enough (during 1 billion years) formation of such sys- 
tems a.t z ^ 4 (Pipino & Mantcncci 2004), which 
enables application of photometric methods to their 
investigation. The efficiency of selecting such galaxies 
with the aid of radio astronomy methods, beginning 
from moderate redshifts {z > 0.5) is confirmed by 
several groups (Pedani 2003). A combined diagram 
of Hubble "K-z" for radio galaxies and field galax- 
ies (Jarvis et al. 2001; De Breauk et al. 2002) show 
that radio galaxies have the highest luminosity at any 
redshift < z < 5.2 (Reuland et al., 2003). Besides, 
radio galaxies have supermassive black holes whose 
mass is generally proportional to a stellar bulge one 
{Mbh ~ 0.006Mb„;d<,e, Maggorrian et al. 1998), and 
this fact is additional evidence of the presence of an 
already formed stellar population. Note that the esti- 
mate of the age of distant galaxies is also of interest 
in connection with searching for primeval black holes 
with masses 10^ - 10^ Mq. 

Formation of radio galaxies at redshifts z ^ 3 — 5 
provide the already formed stellar populations at 
~ 2 — 4 in the ACDM models. Thus, when selecting 
distant radio galaxies, wc isolate with sufficient effi- 
ciency giant elliptical galaxies which can be used to 
estimate the age of a stellar population and to inves- 
tigate the ratio t{z) (Parijskij 2001; Verkhodanov & 
Parijskij 2003, Starobinsky et al. 2004). 

The present paper describes approaches, methods 
and results concerning cosmological parameters es- 
timates using samples of elliptical galaxies. We will 
discuss, first, the problems connected with the use of 
evolutionary synthetic models of spectra of galaxies. 
Then, using the data on the "Big Trio" project (Par- 
ijskij et al. 1996, 2000a, 2000b; Parijskij 2001; Sobol- 
eva et al. 2000; Kopylov et al. 1995; Verkhodanov 
et al. 2002) and other authors (Verkhodanov et al. 



1999) we will present results of the first attempts of 
correction of the standard cosmological model with 
cold dark matter (SCDM) by means of age estimates 
of parent galaxies responsible for the origin of pow- 
erful radio galaxies at large redshifts. Further we will 
give a summary of attempts to estimate the relation- 
ship between the age of galaxies and their redshifts 
from the current evolutionary models of the stellar 
population of elliptical galaxies for a wider interval of 
redshifts, including close, z < 1, objects. 

2. Photometric dating 

2.1. Evolutionary models of spectra 

At the end of the 1980s and in the early 1990s at- 
tempts were made to employ color characteristics 
of radio galaxies to estimate redshifts and ages of 
stellar systems of parent galaxies. Numerous evolu- 
tionary models were proposed, which led to result 
strongly different from one another in comparing with 
observational data (Arimoto & Yoshii 1987; Cham- 
bers & Chariot 1990; Lilly 1987,1990; Parijskij et 
al. 1996). Over the past few years the models PE- 
GASE: Project de'Etude des Galaxies par Synthese 
Evolutive (Fioc & Rocca-Volmerange, 1997) and GIS- 
SEL'98: Galaxy Isochrone Syntheses Spectral Evolu- 
tion Library (Bruzual & Chariot 1993; Bolzonella et 
al., 2000), in which the defects of previous models 
were eliminated, have been widely used. 

In the "Big Trio" experiment (Parijskij et al. 1996) 
we also attempted to apply these methods to dis- 
tant objects of the RC catalog with steep spectra, 
for which we measured the values in the four filters 
(BVRI). The procedure of smoothing was used, which 
made it possible to simiilatc and predict the flux in 
the given filter with the given SED with allowance 
made for the filter response function of this band and 
also with the effects of the redshift allowed for. These 
changes in the procedure permitted the results to be 
more reliable in comparison with the previous paper 
work. 

Preliminarily (Verkhodanov et al. 1999) we inves- 
tigated applicability of new models to populations of 
distant (z > 1) radio galaxies with the known red- 
shifts, for which we have managed to find in the lit- 
erature more or less reliable data of multicolor pho- 
tometry in the optical and near-infrared not less than 
in three filters. In particular, it is shown that red- 
shifts can be estimated to an accuracy 25 — 30% at 
1 < z < 4, given the measured stellar magnitudes in 
more than three filters. But if, at least, one bright- 
ness estimate in the infrared range is available, than 
it suffices to use measurements in the tree filters. Es- 
timations were made for two evolutionary models PE- 
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GASE (Fioc & Rocca-Volmerange 1997), which was 
constructed for the galaxies of the Hubble sequence 
both with star formation and passively evolving. One 
of the advantages of this model consists in the exten- 
sion to the near-IR (NIR) range of the atlas of the 
synthetic spectra of Rocca-Volmerange and Guider- 
doni 1988). This model reconsiders a library of stellar 
spectra which is computed with allowance made for 
parameters of cold stars. The model covers a range 
from 220A to 5 microns. According to the authors, 
the algorithm of the model traces rapid evolutionary 
phases, such as those of the red supergiant or AGB 
in the near-IR range. For the computation a wide set 
of SED curves was used for massive elliptical galaxies 
in a range of ages 7 x 10^ years to 19 x 10^ years. 
We have also used the computations for the ellipti- 
cal galaxies of the library of synthetic spectra of the 
model GISSEL'98 (Bolzonella et al. 2000). The spec- 
tra are constructed with the aid of the evolutionary 
models of Bruzual and Chariot (1993,1996). For the 
calculation of the synthetic spectra of the elliptical 
galaxies of this library, the following parameters of 
star formation were assigned: simple star formation 
(SSP - simple stellar population) , the duration of the 
process of star formation is 1 billion years, while de- 
caying of the burst activity of star formation proceeds 
by an exponential law. The model used solar metal- 
licity. The initial mass function (IMF) with an upper 
limit of 125 solar masses has been taken from the pa- 
per by Miller and Scalo (1979). As is shown in the 
paper by Bolzonella et al. (2000), the choice of the 
IMF does not effect the accuracy of determination of 
redshifts. The model tracks are calculated in a wave- 
length range from 200 to 95800A. For our computa- 
tions we have used the range assigned by a redshift 
limit from to 6. The sets of evolutionary models are 
accessible at the server ,http : //sed. sao . ru (Verkho- 
danov et al., 2000). 

2.2. Procedure of estimating the age and red- 
shift 

Prior to the application of model curves we carried 
out their smoothing by the filters with the application 
of the following algorithm (Verkhodanov et al., 2002): 



n/2+jJj 



(3) 



where Si is the initial model SED curve, Sik is the 
one smoothed by the fc-th, /fc(z) is the curve of trans- 
mission of the fc-th filter "compressed" (1 -I- z) times 
when moving along the axies of the point in the fil- 
ter response function. From the k curves of SED thus 



formed, we have constructed a two-dimension array 
(A-filter) of smoothed synthetic stellar spectra for fur- 
ther computations. 

The estimation of ages and redshifts of radio 
galaxies was made by the method choosing on the 
smoothing curves of SED of optimum positions of 
photometric values obtained in different bands in the 
observations of galaxies. We have used SED curves 
already computed and stored in tables for different 
ages. The algorithm of the choice of optimum posi- 
tion of points on the curve consisted (Verkhodanov 
1996) in shifting the observational points along the 
SED curves. In so doing, such a position was found 
at which the sum of the squares of the departures of 
the points from the corresponding smoothed curves is 
minimum, i.e. the minimum of was actually com- 
puted 



N filters 

E 

fe=i 



Fobs,k -p-SEDk{z) 



(4) 



where -Fobs , k is the observed stellar magnitude in the 
fc-th filter, SEDfc(z) is the model stellar magnitude for 
the given spectral distribution in the fc-th filter at the 
given 2, p is the free coefficient, ak is the measurement 
error. The redshift was determined from the shift of 
the position of the observed magnitudes at their best 
position on the SED curves from the position "rest 
frame" . From the general set of curves we chose such 
ones on which the sum of the squares of discrepancies 
for the given observations of radio galaxies prove to 
be minimum. 

We checked the correctness of estimates of ages 
(and redshifts) by 2 methods. In the first one we took 
synthetic spectra obtained by means of smoothing by 
filters the SED curves for different age. This proce- 
dure made it possible to simulate CCD observations 
for 5 filters. Further the points were chosen corre- 
sponding to the filters VIJHR, for instance, at the 
redshift z — 0.54 and also the model GISSEL with 
SEDs for 1015.1 and 5000 Myr. Two tests were ap- 
plied for each age to estimate the value: with fixed 
z — 0.54 and unfixed redshift. From the results of 
testing a conclusion can be drawn that the age and 
redshift are defined reliably, however, a falling on 
neighboring curve of ages is possible, which gives an 
error of 200 Myr, while at unfixed z the result is also 
affected the discretization in wavelength A in the SED 
curves (the error in z is up to 60%). 

In the second case capabilities were studied of the 
method of determination of the redshifts and ages 
of the stellar population of parent galaxies from the 
data of multicolor photometry. For this purpose, we 
have selected about 40 distant galaxies with known 
redshifts, for which stellar magnitudes in not fewer 
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than four filters (Verkhodanov et al., 1998b, 1999) 
are presented in the hterature. At first, using the se- 
lected photometric data with the use of the models 
PEGASE and GISSEL'98 only ages of the stellar pop- 
ulation of parent galaxies at a fixed known redshift 
were determined. Then a search was made for an opti- 
mum model of the SED curve with a simultaneous de- 
termination of the redshift and age of the stellar pop- 
ulation. After that a comparison of the obtained val- 
ues was made. By this method we estimated both the 
age of the galaxy and the redshift within the frame- 
work of the given models (see also Verkhodanov et al. 
1998a, 1999). It is clear from general considerations 
that the reliability of the result at large redshifts is 
strongly affected by the presence of infrared data (up 
to the K range) because when fitting we overlap the 
region of rapid jump of the spectrum before the region 
of SED, and thereby we can reliably, with a well de- 
fined maximum on the plausibility curve, determine 
the position of our data. Indeed, when checking the 
reliability of the procedure with the aid of the mea- 
surements available with keeping of only 3 points, one 
of which is in the K range, we obtain the same result 
on the plausibility function as from 4 or 5 points. If 
the infrared range is not used, then the result turns 
out to be more uncertain. However, as it is shown in 
the paper by Verkhodanov et al. (1999), the variant 
of 4 filters close disposition as in our case of BVRI 
photometry yielded a good result in the sample of 6 
objects. This result coincides with the one obtained 
when all the filters were used, including the infrared 
range. 

3. Sample of objects 

It should be noted that the choice of elliptical galax- 
ies as objects for our investigation is not accidental. 
They can be considered as the most optimum objects 
among stellar systems having a homogeneous enough 
stellar population. Although such objects have (mod- 
erate) metallicity gradients (Friaca & Terlcvich 1998), 
the modeling showed (Jimenez, Loeb, 2002) that the 
variation of metallicity leads to an uncertainty of es- 
timates of the age within 0.1 Gyr, which lies inside 
uncertainties of estimates. 

In the given investigation we use radio galaxies 
which, as a rule, are identified with giant elliptical 
galaxies (gE) and are good "lanterns" and representa- 
tives of distant stellar systems. The standard point of 
view of the last decades has been that powerful radio 
galaxies are associated with old huge stellar system 
of the gE-type having the red color. The experience 
of using globular clusters in our Galaxy to estimate 
the age of the Universe shows that the search for the 
oldest stellar systems at large redshifts may be useful 



for chronometrying of the rates of expansion of the 
Universe at any distances at which powerful radio 
galaxies still existed. As many groups (Rawlings et 
al. 1996; van Beugel et al 1999), including the project 
"Big Trio" (Soboleva et al. 2000) have shown that 
powerful galaxies appeared at redshifts of about 5. 
The whole interval < z < 5 can be potentially in- 
vestigated even today since the sensitivity of radio 
and optical telescopes is sufficient for investigation 
of such powerful radio and optical objects. In con- 
trast to quasars, the radiation of the stellar popula- 
tion can be readily isolated from that of the gaseous 
component. Note, however, that in radio galaxies un- 
certainties arise in photometric determination of the 
age because of different factors (see, for example, Moy 
& Rocca-Volmerange, 2002), such as ionization and 
transillumination of the radiation from the nucleus, 
interaction of clouds and jets etc. Besides, the galax- 
ies at early stages could be interacting, which changes 
the stellar population. Nevertheless, the radio galax- 
ies remain so far the only simple means of studying 
elliptical galaxies at large redshifts. 

3.1. Data on radio galaxies from the catalog 

"Cold" 

The given sample is composed from FRII-typc galax- 
ies found in the RATAN-600 survey "Cold" (Pari- 
jskij et al.. 1991, 1992) with involvement of data of 
multicolor photometry for estimating color redshifts 
and ages of stellar systems of parent galaxies (Par- 
ijskij et al. 1996; Verkhodanov et al. 2002). Later, 
spectral observations at BTA with the device SCOR- 
PIO (Afanasiev et al. 2002) were carried out, which 
confirmed with high accuracy (correlation coefficient 
0.92) the photometric estimates. 

In the program "Big Trio" BVRI values of about 
60 radio galaxies were estimated, and it was discov- 
ered that although the color age docs not hav(^ a large 
dispersion, the upper limit of the age is a sufficiently 
reliable function of redshift (the larger z, the less the 
maximiim age). A comparison of this upper limit with 
the SCDM model showed that age is not at variance 
with the SCDM model without the cosmological con- 
stant A, but in the interval 0.7 < z < 2 there are 
objects with the color age greater than the age of the 
Universe at the corresponding redshift. Such a sit- 
uation, as it is known (see, for instance, Sahni and 
Starobinsky 2000) is eliminated in the ACDM mod- 
els. Indeed, the age of such a Universe does not differ 
from the SCDM model either at very small or at very 
large redshifts, which is seen from the formulae pre- 
sented in the paper mentioned. 

However, in the interval of redshifts z = 1 — 2 the 
difference may reach 1-2 billion years, which is close 
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Figure 1: Histogram of the distribution of galaxies in z with a formal age above that of the Universe. 



to the possibilities of the experiment. The simple the- 
ory A 7^ for a spatially flat isotopic cosmological 
model yields a relationship between the position of 
a maximum of departures from the SCDM model on 
the axis of redshifts and the value of cosmological con- 
stant. The first attempts to estimate the value of the 
cosmological constant from the "Big Trio" data were 
made in 1999 (Parijskij 2001). In connection with the 
uncertainty in quantitative estimates of the age of 
galaxies from the measured redshifts, percentage of 
galaxies whose age formally exceeds that of the Uni- 
verse in the SCDM model with A = was estimated. 
Then a histogram of distribution of number of these 
galaxies as a function of z was constructed (Fig. 1). 

From the position of the maximum the proportion 
of "dark matter" (Ha) which turned out to be close 
to the value derived from la-type supenovae was esti- 
mated: 57a — 0.8 — 0.6 (Fig. 2). This result stimulated 
our further steps in the usage of age characteristics of 
the stellar population. Fig. 3 shows all the data col- 
lected by the group "Big Trio" in 2001 from steep 
spectrum radio galaxies. It is seen that there present 
objects with large z, but at z < 2 dispersion of ages 
is great. But the larger the redshift, the less the age 
of the oldest object, as it was to be expected in all 
evolutionary models of the Universe. Having chosen 
a population of objects the age of which is younger 
than that of the Universe in the ACDM model by 
more than 2 billion years, obtain a relationship z(t) 



similar to the one displayed in Fig. 4 The relation- 
ship of such a type can already serve as a basis for 
estimating R(t). 



3.2. Data on investigated radio galaxies with 

z > 1 

As it was stated above, to check the procedure and 
estimate redshifts and age of stellar systems, we made 
a sample of radio galaxies FRII with redshifts up to 
z = 3.80 (Verkhodanov et al. 1998b, 1999) using data 
obtained by other authors. 

It should be noted that literature photometric 
data are very inhomogeneous. They were obtained not 
only by different authors but also with the help of dif- 
ferent instruments with different filters. It was not al- 
ways that measurements for one and the same object 
were made with equal apertures etc. For this reason, 
after the final selection from 300 radio galaxies of the 
initial sample only 42 remained. The greater part of 
the sample turned out to be beyond the limits of the 
sample because they have the properties of quasars, 
which impedes strongly the use of the procedure SED 
for standard elliptical galaxies The derived median 
value of the age for the given sample is 5 Gyr for the 
model GISSEL and 9 Gyr for the model PEGASE. 
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Figure 2: Estimates ofil\ from the objects of the project "Big Trio' 
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Figure 3: Relationship t{z) for radio galaxies having steep radio spectra with large z taken from different 
published papers. It is seen that at z <2 the dispersion of ages is large. 



3.3. Clusters of galaxies 

The subsample of elliptical galaxies from clusters, 
which was proposed by A. Kopylov (2001) is the most 
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Figure 4: Relationship z{t) of objects with tstarform — tuniver < 2 billion years 



representative from the investigated group of objects 
(Table 1). We have used for its compilation the data 
from the paper by Stanford ct al. (2002) contain- 
ing a sample of 45 clusters of galaxies at redshifts 
0.1 < z < 1.3. For all the objects photometric data 
from the optical and ncar-infrarcd region arc avail- 
able. On the average, stellar magnitudes in the bands 
BIJHK are presented for each galaxy. For our sam- 
ple wc have selected by the color index 5-9 objects, 
typical elliptical galaxies, from 25 clusters, a total of 
175 elliptical galaxies. Table 1 presents the selected 
clusters with the numbers of galaxies which photo- 
metric data are used to estimate the age of stellar 
systems. The bands in which observations of clusters 
were carried out, their redshifts and K-values are also 
presented. 

4. Procedure of estimating parameters 

Our approach is based on the analysis of the function 
t{z): 

CO 

z 

constructed from ages of radio galaxies depending on 
the redshift. As the function H{z) we used the ex- 
pression 

= Hl[rL^{l + zf + A + B{1 + z) + C(l + zf], (6) 



where A + B + C = 1- 

The fitting of function t{z) to the data analyzed 
was performed with the aid of variation of four pa- 
rameters {Hq, f^TO! A, B). We divided the whole set of 
redshifts into equal intervals lS.z and used the max- 
imum age value in each of the intervals. From the 
sum of the squares of discrepancies a four-parameter 
plausibility was constructed. With the values of the 
parameters _B C = 0, i.e. when the simplified model 
of the function H{z) defined only by two parameters 
{Hq, Q.m) and A = 1 — fl^ = Oa was used, there 
arc stable solution of both models of evolution of the 
stellar population. The results of determination of the 
parameters are listed in Table 2, in which the param- 
eters of approximation of the curve for the intervals 

= 0.2 and 0.3 are given for also in Fig .5. 

4.1. Effects the errors have on the estimates of 
peirameters 

This method of determination of Hq and is stable 
enough concerning the input parameters and system- 
atic effect. As the modeling has shown the variation 
of the initial metallicity caused a change in the age 
by 0.1 Gyr (Jimenez & Loeb 2002). The change in the 
initial mass function does not affect the model SED 
either (Bolzonella et al. 2000). 

The error in the determination of the age, which 
may be connected with wrong classification of the 
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Tabic 1: Selected elliptical galaxies members of clusters. The table contains the following fi,elds: equatorial 
coordinates and name of the cluster, numbers of selected galaxies from the cluster according to the published 
order (Stanford et al. 2002), redshift, used filters, K-limit 



KA+Uec(2000) 


Name 


ID numbers of galaxies 


z 


niters 


TT- 


001631.2+791649 


3C 6.1 


14,18,33,34,43 


0.840 


KJIR 


19.8 


001833.5+162515 


CI 0016+16 


8,9,15,17,22,27,35,40,42 


0.545 


KHJIV 


19.1 


002354.5+042313 


GHO 0021+0406 


9,12,15,24,29 


0.832 


HJIR 


20.0 


002635.7+170945 


/^l AAO A \ 'X f 

CI 0024+16 


3,4,5,8,9,10,16,26,33 


0.391 


KHJRg 


18.8 


UU49 1 U . 9 — z44U43 


Vidal 14 


6,11,14,17,23,31,40,43 


A C O A 

0.520 


KJiV 


1 O A 

18.0 


005657.1—274030 


T1 OOO 1 C/~^T 

J1888.1dCL 


7,17,23,45 


A CT/? A 

0.560 


KHJiV 


19.2 


011018.5+314719 


3C 34 


7,8,13,21,23,34,40 


0.689 


KHJiV 


19.1 


030619.1+171849 


CHO 0303+1706 


6,7,13,18,22,24,34,40 


0.418 


KHJRg 


18.8 


A O O A A 1 /I 1 1 COOAA 

032001.4+153200 


CrHU 0317+1521 


o o r\ io 1/1 OO 0/1 

3,8,9,13,14,23,24 


A COO 

0.583 


T/" TTT T 

KJIV 


19.2 


041246.6—655055 


1 1557. 191C 


OA or OA OT OA 

20,25,30,37,39 


0.510 


KHJIV 


19.1 


A /I C ^1 1 A A A OA A C T 

045410.9—030057 


MS 0451.6-0306 


Oil 10 1AOC01 00/1A 

3,11,12,19,25,31,33,40 


A CT O A 

0.539 


KJlJiV 


1 A O 

19. 2 


A^O A A /I O 1 TAOOl C 
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KJI 
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093239.6+790632 


30 220.1 


5,8,12,16,17,19,24,25 


0.620 


KHJIV 


19.5 


105659.5-033736 


MS 1054.5-032 


4,7,9,13,14,22,23,25,26,30 


0.828 


KHJiR 


20.3 


114022.2+660813 


MS 1137.5+6625 


3,8,10,11,12,15,16,17,20,21,24 


0.782 


KHJiR 


20.0 


132448.9+301138 


GHO 1322+3027 


2,5,8,10,17,18,20 


0.751 


KHJiR 


20.3 


141120.5+521210 


30 295 


9,18,24,27,28,31,32,34,35,37 


0.461 


KHJiV 


18.8 


151100.0+075150 


30 313 


6,12,13,21,33 


0.461 


KJiV 


18.5 


160312.2+424525 


GHO 1601+4253 


4,6,17,18,19,22,43 


0.539 


KHJiV 


19.2 


160424.5+430440 


GHO 1603+4313 


5,15,17,36,40 


0.895 


KHJiR 


20.3 


160436.0+432106 


GHO 1604+4329 


7,10,20,21,28,30,32 


0.920 


KHJiR 


20.1 


205621.2-043753 


MS 2053.7-0449 


35,39,51,89 


0.582 


KJIV 


19.2 


220403.9+031248 


GHO 2201+0258 


10,11,13,14,17,25,29,35,38 


0.640 


KJIV 


19.3 



Tabic 2: A two-parameter fitting of cosm,ological parameters by formulae (5) and (6) (at B = C = O.Oj for the 
approximation curves in the interval Az — 0.2 and 0.3 for both m.odels of stellar population. In the columns 
are presented: the used model of the stellar population, the interval flm, ^A, Ho, the discrepancy e of the 
relationship t{z), the relative discrepancy c/Tq, where Tq=13.7 Gyr is the age of the Universe. 



model 






Oa 


Ho 


e 




SED 










[Myr] 




GISSEL 


0.2 


0.2 


0.8 


77.7 


1695 


0.12 


GISSEL 


0.3 


0.2 


0.8 


71.5 


1367 


0.10 


PEGASE 


0.2 


0.2 


0.8 


65.4 


4101 


0.30 


PEGASE 


0.3 


0.2 


0.8 


53.0 


2748 


0.20 



galaxy type and, therefore, with the choice of SED, 
plays in this case in favor of the approach being de- 
scribed. The models SED corresponding to elliptical 
galaxies give the oldest stellar population among all 
the spectra. Thus, if in the galaxy are found star for- 
mation regions giving a contribution to photometric 
data, then the choice of the optimum model is dis- 
placed toward a younger stellar population, and the 
galaxy falls out from our sample. 



To check the influence of the sample completeness 
on the result, we used the bootstrap method consist- 
ing in multiplying of the initial sample and in con- 
structing a new one by means of accidental choice 
of objects. The coeflicient of multiplying the sample 
was taken to be cqiial to 100, and for each interval 
IS.Z the number of objects equal to the initial num- 
ber was chosen in an accidental manner. 50 tests were 
thus applied, and in each case the values of the pa- 
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PEGASE model (dz=0.2): 

4 parameters fitting 

2 parameters fitting 




PEGASE modei (dz=0.3): 

4 parameters fitting 

2 parameters fitting 






GISSEL model (dz=0.2): 




4 parameters fitting 




2 parameters fitting 















GISSEL model (dz=0.3; 



4 parameters litling 
2 parameters fitting 




Figure 5: The relationship t{z) for models given in Table 2. Upper figures estimation for the model PEGASE, 
low figures — the models GISSEL. Left figures — discretization of Az=0.2, right ones — discretization of 
Az=0.3. The curves are calculated for parameterization (5,6) using maxima of galaxy ages in the given redshift 
intervals. 



rameters were estimated. As a result, dispersion of 
the estimates = 72 ± 7 for the model GISSEL 

and iJo = 53 ± 6 for the model PEGASE for the in- 
terval Az = 0.5 were obtained. We should note that 
accuracy of ±7 for the Hubble parameter is the inner 
accuracy of the method, and actually one elder galaxy 
than existing in the sample can shift the estimate, but 
not far. 

The values of remained unchanged, which is 
explained by the influence of objects at moderate red- 
shifts z = 0.3—1.0 having a relatively small dispersion 
of ages. We disregarded the effect of the interval self- 
absorption in the galaxy (Sokolov et al. 2001) because 
of the ambiguity of solutions for a small number of 
the used input parameters (number of filters) in our 
case and the necessity for the determination of a large 
number of unknown parameters. 

The possible great contribution to the discrepancy 
in the approximation of the relationships, which pro- 
diicc^s the effect of accidental errors, can be consid- 
erably reduced by further statistical accumulation of 
data since the number of known galaxies of type FRII 
will amount by different estimates to a few thousand 
in the near future. 



5. Discussion of results 

The results of our paper confirm, firstly, that we live 
in an evolving Universe. Secondly, that in contrast 
to the standard relativistic flat model AGDM is situ- 
ated within the errors of the method. To estimate the 
quintessence ^(^(z), the accuracy is not sufficient yet. 
For the united data of different populations of ellip- 
tical galaxies, for radio galaxies too, an analysis was 
made of the upper limit of the age of formation of stel- 
lar systems. Prom these data boundaries of determi- 
nation of the cosmological parameters Hq and A term 
were estimated: i?o = 53 ± 10, and f^A = 0.6 ±0.1 
in the model PEGASE. Note that the models GIS- 
SEL have lower dispersion of ages for each interval, 
i.e. they give a more stable result and, possibly, as a 
consequence, more reliable cosmological parameters. 
The discrepancies of the relationship t{z) decrease 
and, therefore, the accuracy of estimates of param- 
eters, when changing from the interval Az = 0.2 to 
the interval Az = 0.3, improves, which is explained 
by more reliable determination of the maximum age 
at a larger interval. 

As far as the procedure applied is concerned, it 
should be noted that one of the main problems is 
the use of investigations of radio galaxies in which. 
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apart from standard evolution of stars, the photomet- 
ric measurements may be affected by other factors 
as well. Nevertheless, new models (for instance, PE- 
GASE2: Le Borgue & Rocca-Volmerange 2002) which 
allow these factors to be taken into account, are al- 
ready beginning to appear. 
Note also that in principle 

1) use of age characteristics of galaxies for in- 
dependent estimates of cosmological parameters is 
prospective; example — the estimate of the A-term, 
which can be improved by extending the sample and 
employing more refined models; 

2) the first examinations of the test "Age of radio 
galaxies — redshift" yielded estimates close to the 
most accurate measurements of parameters from the 
WMAP satellite data (Spergel et al., 2003). 

3) preliminary selection of candidates must be 
performed by different ways since none of the known 
ones is ideal. The relatively powerful radiation (the 
ratio of radio luminosity to optical) points to the fact 
that we deal with a giant galaxy with a supermas- 
sive black hole in the center, which requires time to 
be formed. Unfortunately, there is no generally ac- 
cepted theory of their formation. The proposal to use 
objects in well understood clusters also seems to be 
attractive (Kopylov 2001). Because it is exactly from 
clusters that data for R{t) at small redshifts were ob- 
tained. The use of the upper values of the age in the 
samples at different redshifts seems to be justified. 
Even one the oldest object in the sample is decisive 
in the determination of the lower limit of the Universe 
age at the given redshift, similar to that as the only 
old star in the Galaxy (or the oldest globular cluster) 
defines the minimum age of the Universe today. In 
our opinion the most involved problem is the prob- 
lem of the theory of evolution of the distribution of 
energy in the spectrum of galaxies, and here we have 
considerably various readings. 
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